function y = funcKxfull(nx, nz, zeta, CONSTS)

    k0 = CONSTS.k0;
    c = CONSTS.c;
    Z0 = 4*pi/c;
    
    nxy = func_nxy_from_nz(nz, CONSTS);
    [J1 J2 J1_1 J2_1] = func_J_k(nxy, nx, CONSTS);
    [B1 B2] = func_B_k(nz, nx, CONSTS);
    delta = func_delta_from_nz(nz, nx, CONSTS);
    
    y = ((1i*Z0*k0)/(2*pi))*(J1.*B1./delta + J2.*B2./delta).*exp(-1i*k0*nz.*abs(zeta));
    %y = 2*((1i*Z0*k0)/(2*pi))*(J1.*B1./delta + J2.*B2./delta).*cos(k0*nz.*abs(zeta));
        
end